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Abstract 

Comparison results are given for time-inhomogeneous Markov processes with respect to function classes 
induced stochastic orderings. The main result states comparison of two processes, provided that the com¬ 
parability of their infinitesimal generators as well as an invariance property of one process is assumed. The 
corresponding proof is based on a representation result for the solutions of inhomogeneous evolution problems 
in Banach spaces, which extends previously known results from the literature. Based on this representation, 
an ordering result for Markov processes induced by bounded and unbounded function classes is established. 
We give various applications to time-inhomogeneous diffusions, to processes with independent increments 
and to Levy driven diffusion processes. 
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1 Introduction 

Stochastic ordering and comparison results for stochastic models are topics which have undergone an inten¬ 
sive development in various areas of probability and statistics such as decision theory, financial economics, 
insurance mathematics, risk management, queueing theory and many others. There are several approaches 
for comparing homogeneous Markov processes in the literature. An approach which investigates the in¬ 
finitesimal generator of Markov processes in order to derive comparison results was established by Massey 
(1987). A diffusion equation approach is given for the study of stochastic monotonicity in Herbst and Pitt 
(1991). Bassan and Scarsini (1991) consider partial orderings for stochastic processes induced by expecta¬ 
tions of convex or increasing convex (concave or increasing concave) functionals. For bounded generators 
and in the case of discrete state spaces Daduna and Szekli (2006) give a comparison result for the stochastic 
ordering of Markov processes in terms of their generators. Riischendorf (2008) established a comparison re¬ 
sult for homogeneous Markov processes using boundedness conditions on the order defining function classes. 
Comparison results for homogeneous Markov processes with transition functions defined on general Banach 
spaces are given in Riischendorf and Wolf (2011). The results are based on an integral representation of 
solutions to the inhomogeneous Cauchy problem. The directionally convex ordering of a special system of 
time-inhomogeneous interacting diffusions was considered in a similar way in Cox et al. (1996) and Greven 
et al. (2002). 

Mainly motivated by financial applications a stochastic analysis approach has been developed in El Karoui 
et al. (1998), Bellamy and Jeanblanc (2000), Gushchin and Mordecki (2002) and Bergenthum and Riischen- 
dorf (2006, 2007a). In these papers comparison results for d-dimensional semimartingales are established 
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based on the Ito-formula and on the Kolmogorov backward equation (see also Guendouzi (2009)). A coupling 
approach for diffusion processes and stochastic volatility models has been developed in Hobson (1998). An 
approximation method is used in Bergenthum and Riischendorf (2007b) to give some comparison results for 
Levy processes and processes with independent increments. Several examples and applications to ct-stable 
processes, NIG processes and GH processes are discussed. For multidimensional Levy processes using an 
analytical formula Bauerle et al. (2008) investigate dependence properties and establish some comparison 
results for the supermodular order. They study the question, whether dependence properties and orderings 
of the distributions of a Levy process can be characterized by corresponding properties of the Levy copula. 

Comparison results for time-inhomogeneous Markov processes based on the theory of evolution systems 
on general Banach spaces as used in this paper have not been investigated before. For our main comparison 
result we establish a representation result for solutions of the evolution problem associated with a family of 
infinitesimal generators. We do not use for our approach approximation arguments or coupling arguments as 
in Hobson (1998), Greven et al. (2002) or in Bergenthum and Riischendorf (2007b). Moreover, the application 
of the theory of evolution systems on Banach spaces allows us to reduce regularity assumptions necessary in 
the stochastic analysis approach based on Ito’s formula. 

Applications of this comparison result are given to processes with independent increments (in the sequel 
abbreviated as PII), inhomogeneous diffusions and to diffusion models driven by Levy processes. Therefore, 
we introduce generators of several interesting orderings for which the conditions of the comparison result 
do hold. Since we are interested in comparison of Markov processes (X t )t> o and ( Y t )t>o w.r.t. an integral 
stochastic order <jr, that is Ef(X t ) < Ef{Y t ), t > 0 for all / € A, an integrability condition like 

f £ p| C\P Xt ) n C\P Yt ) for all /£J (1.1) 

t> o 

is indispensable and is made from now on. Also for general state spaces E let T be a set of real functions 
on E in some Banach function space B and let <jr denote the corresponding stochastic order on Ad 1 (E, $), 
the set of probability measures on E, defined by 

l-i <j= v if J fdji < J fdu, (1.2) 

for all f £ P such that the integrals exist. Some interesting examples of stochastic orderings <jr are given 
by the following function classes T for E = R d , 


At := {/ 

R d — >• R; / is increasing} 

(1.3) 

Ax := {/ 

R d —5- R; / is convex} 

(1.4) 

-Fdcx := {/ 

R d — > R; / is directionally convex} 

(1.5) 

Am := {/ 

R d —5- R; / is supermodular} 

(1.6) 

T■ — T 

J 1CX • J cx 

Fl At i *^~idcx -Alex Fl -At 5 Asm -Am Fl -Ait ■ 

(1.7) 


or by subclasses of them. Orderings induced by one of these function classes T are also generated by Afl C °°, 
where C°° is the set of all infinitely differentiable functions, as well as by many other order generating function 
classes P° C T. For notions, properties and applications on stochastic orders we refer to Tong (1980), Shaked 
and Shanthikumar (1994) and Muller and Stoyan (2002). 

As mentioned before the overall aim is to compare Markov processes w.r.t. stochastic orderings <jr. 
For stochastic orderings induced by bounded function classes the Banach space B(M. d ) of bounded functions 
respectively C°°(v) of measurable and essentially bounded functions, where v is a suitable measure, respec¬ 
tively Cb(K. d ) is utilized. For processes with translation invariant transition functions comparison results are 
given for unbounded function classes on modified £ p -spaces, which are introduced for this purpose in Section 
2. We establish easy to verify and flexible ordering criteria which allow us to apply these results to general 
classes of models as for example to PH and to Levy driven diffusion processes and to general order defining 
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function classes. Levy driven diffusions, PII, stochastic interest rate models and stochastic volatility models 
driven by a Levy process have found considerable recent attention in the financial mathematics literature. 
Using the general frame of evolution system and Banach function space theory we are able to treat such 
classes of time-inhomogeneous Markov processes. 

In detail our paper is structured as follows: In Section 2 we introduce some relevant notions and results 
from evolution system theory on general Banach spaces. We establish that transition functions associated to 
inhomogeneous Markov processes are evolution systems on certain Banach spaces. In particular we introduce 
£ p -type spaces with modified £ p -norm and prove that the transition operators define evolution systems on 
these modified £ p -spaces. As a consequence several interesting examples like comparison w.r.t. convex 
functions can be dealt with in generality. In Section 3 we introduce the weak evolution problem. The main 
tool of our ordering method is the representation theorem for solutions of the weak evolution problem given in 
Theorem 3.1. This result extends corresponding results previously known in the literature. As consequence 
we obtain a general comparison result for inhomogeneous Markov processes. 

In Section 4 we discuss several applications that can be dealt with by the generalized approach in this 

paper. We apply this approach to Levy driven diffusion on C'o(IBU) in Section 4.1, to PITs in Section 4.2 and 
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to Levy driven diffusions on the unbounded function class C 2 (y) in Section 4.3. In Section 5 we consider 
as particular example the Sobolev Slobodeckii spaces W(M. d ). The transition operators then define pseudo¬ 
differential operators on 'H r (M d ). In this case our representation result is closely related to a representation 
result in Bottcher (2008). 

2 Evolution systems and their infinitesimal generators 

In this section we at first recollect some notions and results from evolution system theory which is strongly 
related to the semigroup theory on general Banach spaces. Our main reference is Friedman (1969); see also 
Pazy (1983) and Engel and Nagel (2000) for examples and applications. 

A two parameter family of bounded linear operators (T Stt ) s <t, s,t € R+ on a Banach space (B, || ■ ||) is 
called an evolution system (ES) if the following three conditions are satisfied: 

(1) T StS = ids = 1, 

(2) T Sit 1 = 1, 

(3) T S ' t = T StU T U)t for 0 < s < u < t. (evolution property ) 

The evolution system is strongly continuous if the map (s,t) T Sjt is strongly continuous for 0 < s < t. T 
is called a contraction, if 

\\T S) tf\\ < ||/|| for 0 < s < t, f Gl. (2.1) 

For a strongly continuous ES we consider the corresponding family of right derivatives or infinitesimal 
generators A s : T>(A S ) C B —> B 


A - f = l sl( T -'+ hf ~ f ) for>> 


(2.2) 


defined on its domain 


V(A S ) := j/ € B lim i ( T StS+h f - /) exists j 


(2.3) 


For the reader’s convenience we give a proof for the following elementary result (cf. Gulisashivili and van 
Casteren (2006), Section 2.3). The corresponding results for semigroups on Banach spaces can be found in 
Dynkin (1965) or Friedman (1969). Integrals on Banach spaces are meant in Riemann sense (cf. e.g. Ethier 
and Kurtz (2005), page 8 and Lemma 1.4). 
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Lemma 2.1 Let (A s , T>(A s )) s> o be the family of infinitesimal generators of a strongly continuous evolution 
system T = (T Sjt ) s < t . Then, it holds: 

(1) For f G B and 0 < s < u < t, 

t-\-h 

lim y [ T s u f du = T s>t f. (2.4) 

/ij,o n J 

t 


(2) If s i—^ T S) tf is right-differentiable for 0 < s < t, then it holds 

d + 

—— T s t f = —A S T S t f. (backward equation) 

ds 

In particular T Stt f G T>(A S ) for 0 < s < t. 

(3) If f G T>(A S ) for s,t € R+, s <t it holds that 

d + 

—r-T s t f = T st A t f. (forward equation) 
dt 


(2.5) 


( 2 . 6 ) 


Proof: Part (1) follows directly from the continuity of t M- T S;t / for all s G [0,f] and / el. For (2) let 
s < t, then the evolution property and the strong continuity of T lead to 


d+ d+ 

-T-T s tf = lim T s s+h —T s t f 
ds Ho ds 

)• rp f rp f T a+h ,tf ~ T a , t f 

= Inn T sa+h —T at f --- 

hi 0 \ds h 

n ,. T s ,tf~ T a<s+h T Stt f 

740 h 


+ lim T ss+h 
/40 


T s +h,tf — Ts,tf 


I\ 


s,s-\-h 


= — lim ■ 

/40 h 

= -A s T s , t f. 


-T a ,tf 


Now let / € T>(Ag), s G [0, £] then we obtain again due to the evolution property 

~dt T = kg h Tet+hl ~ T °.U) 

=Ss r -*(s( T -^-h) 

= T.jAtf 


and part (3) is done. □ 

For a time-inhomogeneous Markov process (X t ) on a measure space (IE,!?) let (P s ,t)s<t denote the 
transition kernel or transition function 

P s ,t(x, V) = P(X t G V | X s = x), x G E, V € S (2.7) 

and T = (T s>t ) s < t the corresponding evolution system, i.e. the transition operator 

T s>t f{x) = Jp St t(x,dy)f(y ) (2.8) 

E 


4 



= E{f{X t ) \X s = x) ( 2 . 9 ) 

for / in a suitable Banach space B of functions on E. This puts Markov processes in the framework of 
evolution systems and evolution equations. 

As first example we consider the Banach spaces C°°{v) := {/ : K. d —> K | / is ^-measurable, \\f\\oo < 00 } 
of z'-measurable and essentially bounded functions, where v is a suitable measure on R d , B( R d ) the class of 
bounded functions on R d and Cb(R d ) the continuous functions in B(M. d ). Recall that for Markov processes 
the properties (1),(2) and (3) from the definition of ES always hold true on these Banach spaces. Moreover, 
observe that uniform continuity of the transition kernels in (2.8) implies strong continuity of the transition 
operator T: 



| P s ,t(x,dy) 


P s ,u(x, dy) | —>• 0 as t —> u, and 



I P s ,t(x,dy) 


Pu,t(x,dy) | —» 0 as s —>• u. 


( 2 . 10 ) 


Proposition 2.2 (ES-property for time-inhomogeneous Markov processes on bounded function 
classes) 

Let ( P s ,t)s<t be the transition kernel of a time-inhomogeneous Markov process (X t )t> 0 on (E, <^). If the 
transition kernel ( P s ,t)s<t satisfies condition (2.10) then the operator defined in (2.8) is a strongly continuous 
contraction ES on B £ {f?(R d ), C°°{v)}. If additionally, f £ C6(R d ) implies T Sjt f £ C^K^), then the 
operator is a strongly continuous contraction ES on B = C&(R rf ). 

Proof: Since T is a family of bounded operators it leaves bounded functions invariant, that is, for / £ Cb(R. d ) 
we obtain 

l|T M /||oo < sup ( Ps,t{x,dy)\f{y)\ < ||/||oo- 

X J 

m d 

Hence, the strong continuity follows due to 

\\T s ,tf — /|| = sup \T S}t f(x) — T StS f(x)\ < 11/11 00 * sup f \P a>t (x, dy) — P s , s (x, dy)\ 

X X J 

R d 


for all s < t and the continuity property (2.10). 

The proof for B £ {_B(R d ), £°°(j/)} is very similar and therefore omitted. □ 


Remark 2.3 (a) For a PII L = (L t )t>o the transition operators on B are given by 

T a , t f(x) = j P Lt ~ La (dy)f (x + y). 

R d 


( 2 . 11 ) 


If B = C b (R d ) the map x T s t f(x) is continuous since the continuity of f transfers to (T s t ) thus 
Cb(R d ) is invariant under (T s j) s <t- In particular the transition kernels of Levy process such as Brow¬ 
nian motion, NIG, VG, GH processes have this invariance property. 

(b) Note that the Brownian motion respectively its associated Brownian kernel 


Ps,t(x, dy) 



R ' 1 


1 

— , =e x P 

\/27r(t - s) 


\x-y \ 2 

2(1 — s) 


dy 
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for x £ K d and s < t is not a strongly continuous ES on Cf,( R“). Hence, it is vital to consider further 
Banach spaces with different norms to establish the ES-property for Markov processes. In Proposition 
2-4 we consider the Banach space of p-integrable functions on R“ and show that in case of translation 
invariant transition function the ES-property does hold. 

(c) In Section 4-1 we will see that Levy driven diffusion defined via the stochastic differential equation (4.1) 
possesses the ES-property on Co- Thus, in order to establish stochastic ordering results induced by 
bounded function classes the Banach space (Co, || • ||oo) is used as a reference space. 

For many stochastic orderings the function spaces Cb(R d ), -B(R d ) or C°°(v) are sufficient; some orderings 
as for example the convex ordering < cx however do not allow bounded generating classes of the ordering. The 
following proposition shows that Markov processes with translation invariant transition functions ( P s ,t)s<t 
are strongly continuous ES on £ p -spaces. Hereto, we consider the p-integrable functions on R d , which we 
denote by 

C p (R d ) = {/ : —> R|/ is measurable, ||/|| p < cxd}, 1 < p < 00. 

Proposition 2.4 (ES-property for time-inhomogeneous translation invariant Markov processes 

on C p (R d )) 

Let (X t )t >0 be a time-inhomogeneous Markov process on with translation invariant transition function 
( P s ,t.)s<t ■ Then the family of transition operators T = (T s ^) s <t defined by 

T s ,tf(x) = J P s ,t{x, dy)f(y) = J P s ,t{ 0, dy)f(y + x) (2.12) 

R d R d 


for f £ C p (M. d ), is a strongly continuous contraction ES on C p (R d ), 1 < p. 
Proof: Let s < t and set P Stt (0,dy) =: P a>t (dy). We have for / £ jC. p (R d ) 


\\T s ,tf — f\\p = 


Ps,t(dy)f{x + y) - /(:r) 


dx 


< J P s ,t(dy) J \f(x + y) — f{x)\ p dx =: J P s>t (dy)h(y), 


where h is a bounded and continuous function with h( 0) = 0 (see Sato (1999, E. 34.10)). For each e > 0 we 
can find a 6 > 0 such that ,<£, P s j(dy)h(y) < e. By the triangle inequality to show strong continuity of 
T s t it is enough to consider for u £ R + the case that s < t and s, t -£ u. Then 


lim [ P s ,t(dy)h(y) < e + lim [ P s t (dy)h(y) 

fu,t—fu, J s—fu,t—fu J 

s <t wd 


= £ ■ 


{|3/I>H 

P u ,u(dy)h(y) 


{M>H 

= £+ f So(dy)hfy) = £. 


{M>H 


Finally, we obtain due to the convexity ofsi-> \x\ p , 


Ts,tf\\l= \ Ps,t(dy)f(x + y) 


dx 


< J P s ,t(dy) J \f{x + y)\ p dx 


(2.13) 
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= J PsMv) J \f( x )\ Pdx = \\f\\p- 

Thus, the operator norm of T Sjt is bounded by one, i.e. ||T Sjt || < 1 and T st f £ £ p (M d ) for all 0 < s < t, i.e. 
T is a strongly continuous contraction evolution system on C p ( R d ). □ 

For some classes of models it is possible to establish the .ES'-property by comparison to the translation 
invariant case. 

Corollary 2.5 Let (X t ) t >o be a time-inhomogeneous Markov process on Ml i with transition function (P s ,t)s<t 
and transition operator T = (T Sjt ) s <t defined by (2.8). If there exists a positive constant c and a translation 
invariant transition function ( Q s ,t)s<t such that 

J f(y)P s ,t( x ,dy) <cj f(y)Q a> t(x,dy), (2.14) 

R d R d 

for all x £ all positive functions f on M d and all s, t £ R + with s < t, then T is a strongly continuous 
ES on £ p (R d ), 1 < p < oo. 

As we see in Proposition 2.4 the translation invariance property of the Lebesgue measure and the trans¬ 
lation invariance of the associated transition function is crucial for T to be a strongly continuous contraction 
ES on £ p (R d ). For a er-finite measure v we circumvent the lack of the invariance property by introducing a 
suitable weighted sup-norm on a sufficiently large subspace of C p (n): 

figM ■= {/ e C p (v) | ||/||p e := sup- f\f{ x + y)\ p dv{x)] P <oo| 

L ^ yGR“ (1+ J/ )* W J J 

R d 

for 1 < p < oo and g > 0. For p = oo we define || ■ ||oo je = || ■ ||oo = || • ||oo,y Thus, the space £^°(^) equals 
The £ p -type space (C p (iy), || • ||p >e ), Q > 0, l<p<ooisa Banach space. 

Lemma 2.6 (C p {v) % || • ||* g ), g>0,l<p<ooisa Banach space. 

Proof: By definition || • ||* is a norm and therefore, £ p Q [y) is a vector space. It remains to show that C p Q {y) 
is complete. This is done similarly to the proof of the Riesz-Fischer Theorem. 

Let (fj)j £n be a Cauchy sequence in C p (u), that is 

II fj - fk\\p >e = sup + || l y ^ e/p ( f I fi( x + V)~ fk( x + y)\ P du(x)j /P 0 (2.15) 

R d 

as j, k approaches infinity. Thus, we have to show that there exists / £ C p e (y) such that fj f for j —> oo 
in C p ( v). Due to 

\\f-fj\\;,g< n/-4ii;, e +ii4-/iii;, e 

it remains to verify this fact for a subsequence of (fj)jeN- We choose the subsequence {fj m )men such that 
( fj m ) converges a.s. and 

OO 

E (2.16) 

m=1 

and denote (/ Jm ) by (f m ) again. Since 

(1 + noii)«/> (/ |9<X+0) '" dv{x) ) s 7(1 ,lr 

R d R d 
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for all g £ C p {v), we have ||/ TO || P < ||/ m ||* e . Due to (2.15) the sequence (/ m ) me N is convergent in C p {v) and 
there exists a limit / £ E p {y). By the Lemma of Fatou we obtain for y £ R d : 

^ + y|) g J I f(x + y)-fm(x + y)\ p dv(x) 

< iimmf ^ + || y ||^ e J I M x + V)~ fm{x + y)\ p dv{x) 

< -fm\\* P ,gY 

l—l oo 

< (5Z ll/fc+t ~ fk\\p, e ) P = (5Z ll/fc+1 - fk\\*p,e) P < 00 ■ 

k—m k—m 

Thus, we have 

S a+wiK- t/ s(Eii/w-aii;,.)- 

By (2.16) the last term tends to zero as rn approaches infinity. This implies the statement of the lemma. □ 

For time-inhomogeneous translation invariant Markov process (^ft)t>o on R d with transition function 
( Ps,t)s<t satisfying a certain integrability condition the ES-property holds for C p Q (v) = C p (v) fl C(K d ), the 
set of continuous functions in C p g (v), with respect to || ■ ||* ,,-norm as well. 

Proposition 2.7 (ES-property for time-inhomogeneous translation invariant Markov processes) 

Assume that ||z|| e is uniformly integrable ui.r.t. P Sjt for some g > 0 with s <t, i.e., 

sup [ \\z\\ s P s , t {dz) —>• 0. 
s<t J K—too 

M >K 

Then, the family of transition operators T = (T s ,t) s <t defined for f £ C p c g (v) by (2.12), is a strongly 
continuous ES on C p cg {v) for 1 < p < oo. If additionally, v is absolute continuous with respect to the 
Lebesgue measure A, then T is a strongly continuous ES on C p g (v). 

Proof: The proof for the strong continuity is analogous to the proof of Proposition 2.4. Hereto, assume v 
to be absolute continuous with respect to A. For / £ C p (u) note that the function 

h*{z) = sup + ^ /p [J \f{x + y + z)~ f(x + y)\ p du(xfj ' P 

is bounded by c(l + || 2 i|| e ) and continuous with h*( 0) = 0 where c = c/ is constant in R + . In the case where 
v is not absolute continuous w.r.t. A, then simply restrict the function class to / £ C p g (v). Consequently, 
h* becomes continuous again. For the boundedness let / £ £ p {v) and observe that the expression 

((i + M)« (1 + b + 2r ) '* 

assumes its maximum in y = 0. Thus, 

h*( z ) < c(sup (1 + || ^ ||)g/p (j I/(* + y + z)\ p dn(x)^j ,P + ||/||;, e ) 



, 1/p 


* ) 

p,e J 


< c ( sup (l±J!i±^! 

_ \ V (1 + ll»ll)* / >' 

■ (l + ||j/ + z||)e/p {j^ + y + *>I’<M*>) 

<c\\f\\;ji + \\z\n 1/p 

for a sufficiently large c > 0. Then we have 

II T s ,tf - f\r p , e = sup (1 + ( J I T s,tf(x + y)- f{x + y)\ p dv{x) S j 

< J P s ,t(dz)h* p (z). 

Note that for each £ > 0 we can find a 5 > 0 such that /{|»|< 5 } Ps,t{dz)h* p (z) < £ due to the continuity 
of h*. Since h* p is bounded by the uniformly integrable function z H > 1 + ||z|| e w.r.t. P Sjt separating the 
integral and letting K.+ , s < t we obtain for M £ K + that 

lim [ P s t .{dz)h* p (z) < e + lim [ P s t {dz)h* p (z) 

" ±u,t—tu, J ’ s—tu, t—yu J 

{\z\>S}r[-M,M\ 

J P s , t (dz)h* p (z) 

{\z\>8}n[—M,M] c 

<£+ [ P u , u (dz)h* p {z) 


s<t 


lim 

s—>u, t—yu 


{|z|><5}n[-M,M] 


lim 

s—>u, t—yu 


P s j{dz)h* p {z) 


{\z\>8}n[-M ,M} C 


<£+ J So(dz)h* p (z) + £ = 2£. 


{\z\>8}n[-M,M] 


Moreover, recall that f Rd \\z\\ e P s ^ t {dz) < oo, thus 

\\T S > t fW™ = S ^P ( i + y| )e /p (/ \ T '*H X + v)\*Mx)) 


i/p 


T(l+IMI) e/p 


(/ J f(x + y + z)P Stt {dz) du(xj^j 


i /p 


R d R d 


< 


T d + IM)«/” (/ P ‘ Adz) ( / 1/(1 + v + z)l " d,/{x) )) 

R d R d 


1/p 


Consequently, we arrive at 

||T s , t /||; e < ||/||* >e sup( j 1 (1 + ||y + «||) e )) 


i/p 


R d 


- I \J lip. 


< C • 


P,,t(dz)(l + \\z\\) e 

R d 

;, e (J Nl e PsAdz)) 1/P <C‘ 


/II f II* 

p,e ‘ 
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Hence, \\T s ,t\c p e (u) II — c ' an< ^ € C p Q (v) for all 0 < s < t, i.e. T is a strongly continuous ES on C p (v). □ 

Remark 2.8 (b— bounded functions) 

The norm modified C p -space, CP(y) allows us to deal with orderings generated by unbounded function classes. 
Let b : R d —>■ [1, oo) be a weight function and define 

B b := {/ : R d -» R | 3c <E R : |/(a;)| < c • b(x)} 

the class of b-bounded functions. In particular, to deal with the convex ordering we choose b{x) = 1 + ||a;|| 
for || • || any norm on R d . Then the class of b-bounded convex functions is a generator of the convex ordering 
< cx and is a subset of C p (u), p < g for any u which integrates (1 + ||a:||) p . 

Thus in the case of C v Q (y) the E.S'-property is implied by the translation invariance property and a 
integrability condition of the transition function ( P s ,t)s<t ■ Under an additional assumption on the transition 
kernel we extend in the following the ES-property to the case of not necessarily translation invariant time- 
inhomogeneous Markov processes for unbounded function classes. More specifically we consider the class 
C\(v) and impose the following assumption ( K ), which is a strengthening of the notion of the Hilbert- 
Schmidt operator: 

We assume that the transition kernel P S}t (s,dz) = k s j(x,x + z)v(dz) corresponding to a time-inhomo¬ 
geneous Markov process has a density w.r.t. v fulfills 

(K) K* t := sup 1 ^ J \k Stt (y,y + z)\ 2 dv(z)j < oo 

R d 


and satisfies the continuity assumption 

(C) For 0 < s < t holds lim K*, t , = K* t . 

For the continuity in (s, s) we make the following local domination assumption. 

(D) For any s > 0 there exists 5 = 5{s ) > 0, C > 0, and a translation in¬ 
variant transition function ( Q s ,t ) such that ||£|| 2 is uniformly integrable 
w.r.t. (Q s ,t) and P s>t f < CQ Stt f, V/ € T\(v), f > 0, £ - s < d. 

Note that conditions (K), (C) are fulfilled under boundedness and continuity assumptions on k. The domi¬ 
nation condition (D) ensuring continuity of in (s, s) might be verified in several applications but could 
also be replaced by further ad hoc assumptions. 


Proposition 2.9 (ES-property for time-inhomogeneous Markov processes on 

Let v be a finite measure with finite second moments on R d and let (X t )t >o be a time inhomogeneous 
Markov process with transition kernel ( P s ,t)s<t which satisfies the assumptions (K), (C) and (D). Then the 
corresponding family of transition operators (T' Sjt ) s < t is a strongly continuous ES on C\(y). 

Proof: For / £ £|( p ) an d € R d holds that 


I T s , t f(x + y)\ 


< 


< 


J k s ,t(x + y, x + y + z)f(x + y + z)v{dz) 

J \k s , t {x + y,x + y + z)\ 2 du(z)^j ( J f 2 {x + y + z)d 


<*) 


1/2 
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This implies that 


We also obtain 


and 


A s ,t ■= sup . 

v 1 + IMI 


-Tttm 


I ( T s , t f(x + y)) 2 dv{x) 

kl t (x + y,x + y +z)dv(z) 
j f 2 {x + y + z)dv(z)J dv(x) 


1/2 


TTbii,/ / 2 (* + » + *><M*> 

1 + ||x + y\\ 1 f £ o / \ 7 /\ 

= 1 + toll i + ll* + «,ll/ /(l + !, + i!) ‘ Mi!) 

<(1+INI)II/Il2,2 


su p x || J kl t (x + y,x + y + z)dv{z) < (1 + ||z||)iC )t . 


As consequence we get 


A,,t< (/(1 + Ikll) 2 ^)) 1/2 ||/Il2 ) 2^f 


Thus by assumption (K) T s t f £ C 2 (v) for / £ G\(y) and by the continuity assumption (C) T Sjt is strongly 
continuous in (s,t) for s <t. By the domination assumption (D) we obtain from Proposition 2.7 that T Sit is 
also continuous in (s, s), 0 < s. □ 


Several modifications of the domination assumption ensuring continuity in (s, s) could be given. 

Example: PII and their infinitesimal generators 

Let L = (L t )t> o be a PII with continuity property (2.10) and with characteristic function given by 

fi e i(t, L t) _ e / 0 ‘ 0s<,i£)ds ^ (2.17) 


where for s > 0 the cumulant function 9 = (0 s ) s >o equals 

+i{£,b s } + J - 1 - z(£,Xc(y)>) F s {dy) (2.18) 

for a cut-off function \c , that is a bounded, measurable real function on R d with compact support and which 
equals the identity in a neighbourhood of zero. Here for each s > 0 the covariance matrix a s is a symmetric, 
positive semi-definite d x d matrix, the drift b s is in R d and F s is a Levy-measure, i.e. a Borel-measure on 
R d which integrates (1 A |cc| 2 ) with F({0}) = 0. 

For a triplet ( b s ,a s ,F s ) defined as above we consider the operator G s f{x ) = G® f(x) + G J S f(x), s £ R+ 
with 


Gs f(x) ■= 7T > F 


E 

j,k =i 


r Jifc. 


d 2 f 

dxjdxk 


(x) 


d 

E 

j =i 


u df, \ 


and 


G J s f(x) 


■= j (f(x + y) — f(x) - 
R d 


d 


E 


-^r{x){Xc{y) 3 )]F s {dy) 


(2.19) 
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on W 2 (V) = {/ G C 2 (R d )n£%(v)\ ^, 9 ^- G £|MV1 < i,j < d} C £|(^). 

For upcoming results it is crucial to know that G s f belongs to C\(v) and whether such an operator can 
be linked to a strongly continuous ES on £ 2 (z/) corresponding to a time-inhomogeneous Markov process. 
Obviously, for / G VV 2 (V) the expression Gf*/ G E\(y) for every s G R+. Moreover, under some mild 
regularity conditions it can be shown that G J S : VV 2 (is) — > F\(y), where W 2 (y) is a subspace of C\{v) with 
norm || • ||$ >2 - 


Lemma 2.10 For the operator G® f introduced in equation (2.19) it holds that Gf*lvv 2 („) : W 2 (^) —> E\(y). 
If additionally 

t 

J J \y\ 2 F s {dy)ds < 00 , 

0 {I y I > 1} 

then the same holds true for the operator G s |^ 2 („). 

Proof: It remains to show that for / G W 2 (v) the jump part G J S f in (2.19) belongs to F\(y). We only 
cover the one-dimensional case. The multivariate version is similar. Let / G VV 2 (^), s G R+, and choose the 
cut-off function Xc(y) = ?yl{|y|<i}- Using the Taylor expansion for x' G R we arrive at 


f(x' + y)—f(x') - f'(x')yl {M< 1 } F s (dy) 


J {fix' + y) - f(x') - f'(x')y)l {M< 1 } F s (dy) 

R 

+ J liV + y) - f(x')\F s (dy) 


{\y\>n 

< J - 2 y 2 \nO\Fs(dy)+ J 

{|y|<i> f|y|>i} 


\f{x' + y)\F s {dy) 


=:h 


— 'I2 


+ \fi x ')\ J (1 A \y\ 2 )F s {dy), 


f|y|>i} 


=’h 


where £ is a suitable intermediate point between x' and x 1 + y depending on y. Thus, the term |/"(£)| is 
bounded on {|y| < 1} since / G G 2 (R d ). For i,zGR, we choose x' = x + z, then, squaring and integrating 
each term successively w.r.t. dv{x). For I 3 we then have 

( f (lA|y| 2 )F 5 (dy)) • f \f(x + z)\ 2 dv(x) 
lM>i} 

< (1 + N)-c 3 sup i+ 1 m (/ \f(x + z)\ 2 dv{x)^j 

<(i + N)-c 3 (ll/li ; >2 ) 2 

for a suitable non-negative constant C 3 . Since 


/ I \ y \ 2Fs ( dy ^ ds < 00 ’ 

0 {|y|>i} 
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for the middle term / 2 we observe that 


f j \f(x + z + y)\ 2 F s (dy)du(x) 

{|i/|>i} 

J (1 + |y+ ^|)( 1 + ^ y + z \ J \f{x + y + z)\ 2 dv(x)jF s (dy) 

{|y+d>i} 

< J (1 + \z + y\) • sup 1 + } l + y] (J \f(x + z + y)\ 2 du{x)^F s {dy) 

<(ll/li;, 2 ) 2 ( J (l + \z\)(l + \y\)F s (dy)) 

{|y|>i} 

<(1 + N)(ll/Il 2 , 2 ) 2 ( / (1 A \y\ 2 )F s (dy) + J \y\ 2 F s (dy )) 

{| 2 /|> 1 } {| 1 /|> 1 } 

<(1 + N)-C 2 (||/Ii;, 2 ) 2 , 

where in the first equality Fubini’s theorem is used, and again c 2 is a suitable non-negative constant. Using 

_2 

that f" £ T 2 (y) we obtain from Taylor expansion that the first term I\ is bounded by a constant ci. In 
summary this yields 


J\G J s f{x + z)\ 2 dv{x) 
R 


J J f(x + z + y)~ /(x + z)~ f'(x + z)yt { \ v \ <1} F s (dy) 


2 

dis(x) 


<(i + H)-(c 1 + C2 (||/||5 i2 ) 2 + c3(||/||5 l2 ) 2 ), 


which implies 


sup 

z 


1 

TTR 


\G J J{x + z)\ 2 dv{x)^j 


< (ci + c 2 


2 ) 2 + C3 


2 , 2 ) 2 ) < 00 . 


Consequently, we proved ||Gf/|| 2 2 < 00 , hence, the statement holds true. 


□ 


Proposition 2.11 (Infinitesimal generator for PII on C\{u)) 

Let 0 be a PII such that E\L t \ 2 < 00 for all t £ M + . Denote by T = (T s ,t) s <t the corresponding 

transition operator on C 2 (v), and its infinitesimal generator A s f defined via equation (2.2) for f £ D(A S ). 
Then, A s f = G s f on VV 2 (^). In particular VV 2 (^) C T>(A S ) for each s £ R+. 

Proof: From Proposition 2.7 we know that T S]t |^ 2 (l/j is a strongly continuous ES on the Banach space 

C\(y). Hence, the limit A s f = lim^o t{T s ,s+hf — f) can be understood in the strong sense on D(A S ). 
Thus, A s : T>(A S ) — > Cl 2 (v). Now, let / £ VV 2 (^) and consider a smooth approximation sequence 1 /at £ 
Cf(M. d ),N £ N such that df]y,d 2 fpf £ C^(R d ) with ||/jv — /|| 2j2 —t 0. Recall that v integrates (1 + ||x|| 2 ), 
thus /at £ jC 2 (u) as well as 3/at, <9 2 /at £ C 2 (v), or equivalently, /at £ V\f 2 (v). From standard literature 2 , e.g. 

1 Utilize standard approximation by a modifier tp which is smooth up to the boundary. 

2 In Dynkin (1965) the infinitesimal generator is given for g E C 2 (M d ) for time-inhomogeneous Levy processes and in Jacob 
(2001) for g E C^(M d ) for Markov processes. 


13 



Dynkin (1965) and Jacob (2001), we know that A s g = G s g for all g G C 2 (R d ), where G s is the operator 
as defined in equation (2.19). Consequently, we obtain A s Jn = G s Jn • The moment condition on (L t )t >o 
yields 

t 

J J \y\ 2 Fs(dy)ds < oo, 
o {M>i> 

Lemma 2.10 then implies 

11^/^115,2 = 11^/^115,2 < \\G?fN\\t 2 +\\G J J N \\l 2 < OO. 

Since G s /at and G s f G L\ {v) for all N G N it follows by majorization 

A s f := lim A s f N = lim G s f N = G s { lirn f N ) = G s f. 

N—>oo N—>oo N—>oo 

In particular, it holds that ||A S /||2 2 < oo and, thus, W 2 (y) C D{A S ) for each s G M+. □ 

Similar results as Lemma 2.10 and Proposition 2.11 hold true under some mild regularity conditions, also 
for time-inhomogeneous pure diffusion processes. 


3 Comparison of Markov processes 

The main goal is to prove a general comparison result for time-inhomogeneous Markov process. Hereto, we 
need a representation result for strongly continuous ES on general Banach spaces. We start by recalling the 
weak evolution problem. 

3.1 Representation result for solutions of the weak evolution problem 

The weak evolution problem is crucial for the comparison result in the next section. Let i be a Banach space. 
For a strongly continuous evolution system T = ( T S}t ) s <t on B with corresponding family of infinitesimal 
generators (A s ) s> o define F t (s) := T Stt f, s <t, f G T>(A a ). Then, by Lemma 2.1 for fixed t, F t is a solution 
of the homogeneous weak evolution equation 

ill ef 

- -Al = -A a u(s) for s < t, (3.1) 

as 

u(t) = /. 

As usual, F t is called the fundamental solution of (3.1). In the autonomous case where A s is independent of 
s, equation (3.1) becomes the homogeneous Cauchy problem. 

Fundamental solutions and the Cauchy problem are studied in Friedman (1969) imposing conditions on 
the resolvent of the generator. In order to deal with inhomogeneous Markov processes we need to study 
solutions of the evolution problem in the following form which is related to Friedman (1969, Part 2, Chapter 
3) (see also Pazy (1983, Chapter 5)). This extension is the foundation of the comparison results in the next 
section. 

For every r G [s,t], s,t G R+ let A r : T>(A r ) C B -> 1 be a linear operator on B and let G(r) be a 
B-valued function on [s, t]. We consider for a B-valued function u on [,s, t], which is right differentiable on 
(s,f), it(r) G D(A r ) for s < r <t, the weak initial value problem or weak evolution problem 

-jp- = -A r u{r) + G(r) for s < r < t, ^ 

u(t) = f. 
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A B- valued function u : [s,<] —> B is a classical solution of (3.2) if u solves (3.2) and is continuous on [s,t]. 

Note that the autonomous case of the weak evolution problem is a weakening of the inhomogeneous 
Cauchy problem. The basic representation result for solutions of the inhomogeneous Cauchy equation in a 
Banach space was used in Riischendorf and Wolf (2011) to establish comparison results for homogeneous 
Markov processes. The following is an extension of this representation result to the weak evolution problem. 

Theorem 3.1 (Representation result) 

Let T = (T S}t ) s < t be a strongly continuous ES on B with corresponding family of infinitesimal generators 
(A s ) s >o- Further, let F t , G : [0, t] —> B for t £ M + be functions such that 

(1) the map r H > T Sjr G(r) is right continuous. 

(2) f* T Str G(r)dr exists for all s,t £ R+ with s <t. 

(3) F t solves the weak evolution problem (3.2), i.e. for s,t £ R + such that s < t, it holds that r F t (r) 
is continuous on [s,t], right differentiable on ( s,t ), F t (r) £ T>(A r ) on (s, t] and 

d+FUs) 

, tK ) =~A 3 F t (s) + G(s) for s<t. (3.3) 


Then, 

t 

F t (s ) = T., t F t (t) - j T Str G(r)dr. (3.4) 

S 

Proof: Let s,t £ R+ with s < t. The right derivatives ^ T s r and ^F t (r) exist due to the assumption on 
T and F t . Using (2.6) and (3.3) we obtain 

L(r., r F,(r)) = VmT.,r+^(F,(r + h) - F,( r)) + Um - T.,,)f,{ r) 

= r,,(^F, W ) + (^T.,)F,W 
= T s , r (-A r F t (r) + G(r)) + T s ,,A r F t (r) 

= T s , r G{r). (3.5) 

Thus (T Str F t (r)) is integrable on [s, t] due to assumption 3.1-(2). Since it is right continuous we can 
integrate (3.5) and obtain 

t t 

J T s , r G(r)dr = J ^-(T s , r F t (r))dv = T a>t F t (t) - T s , s F t (s). 

S S 

Thus, we obtain the representation of the solution of the weak evolution problem 

t 

F t {s) = T a>t F t (t ) - J T s r G{r)dr. □ 

S 

Remark 3.2 (a) Using the theory of pseudo-differential operators, Bottcher (2008) developed an integral 
representation result on certain subspaces of Cq for solutions of an evolution problem related to (3.2). 
See also our Section 5. 

(b) For a special class of infinitesimal generators of strongly continuous semigroups Pazy (1983, Chapter 5, 
Theorem f.2) derived a similar representation formula for evolution systems on Banach spaces. 
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3.2 General comparison result for Markov processes 

We assume that X and Y are two time-inhomogeneous Markov processes with values in (E, S) and let 
S = {S S 't)s<t and T = (T S)t ) s <t denote their strongly continuous evolution systems on some Banach function 
space B on E. Let by A = M s )s>o and B = (B s ) s > 0 denote the corresponding families of infinitesimal 
generator of S and T respectively. Further, we assume that for each s > 0 

Fc V(A s )nV(B s ). (3.6) 

Theorem 3.3 (Conditional comparison result) Assume that 

(1) / £ F implies s H > T Sjt f is right-differentiable for 0 < s < t, 

(2) / £ F implies S s jf £ F for 0 < s < t, (stochastic monotonicity of S) 

(3) A s f < B s f [P' Y °] for all f £ F and s <t. (3-7) 

Then 

S., t f < T„ tt f [P' Yo ] for all f£F,s<t. (3.8) 

Proof: Define for / £ T and t £ R + the function F t : [0, t] —» B by F t (s) := T s>t f — S Si tf ■ Then due to 
Lemma 2.1 satisfies the differential equation 

d+ FAs) 

, y 1 = -B„T a , t f + A s S s t f, F t (t) = 0, on [0,f], 
as 

As consequence we obtain 

d+ 2 (5) = ~Bs{T a ,tf ~ S a , t f) + (A a - B^Ss'tf = -B s F t (s ) + G(s), 

where by (3.6) and assumption 3.3-(2) G(s) := (T s — B s )S s jf is well defined. Further it is non-positive due 
to assumption (3.7). Moreover we obtain that (— B s F t (s )) is also well defined due to assumptions 3.3-(l), 
Lemma 2.1 and 3.3-(2). The continuity of S and T transfers to F t (s). Furthermore, the strong continuity 
of T fl) t and S s>t transfers to the map r i —> T s r G(r), hence, it is right continuous, too. Thus, F t (s) solves the 

weak evolution problem. Since —G(r) > 0 for s < r < t it follows that — /* T s ^G(r)dr exists and is finite 

for all s,t £ M_|_. Thus, the assumptions of Theorem 3.1 are satisfied and imply that the solution F t (s) has 
an integral representation of the form 


t t 

F t (s ) = T ait F t (t) - J T s , r G{r)dr = J T a , r (-G(r))dr > 0, 


as F t (t) = 0 and —G(r) > 0 for all s <t. 


□ 


Remark 3.4 (a) The comparison result for homogeneous Markov processes in Riischendorf and Wolf (2011, 
Theorem 3.1) is implied by Theorem 3.3 since in the case of homogeneous Markov processes the corre¬ 
sponding semigroups fulfil condition 3.3-(l) (see for example Theorem 2.f in Pazy (1983)). 

(b) Let T be any function class in (1.3)-(1.7). If the Markov process X on S. d has a translation invari¬ 
ant transition function (P s ,t)s<t> then X is stochastically monotone ui.r.t. to F(see Lemma 2.8 in 
Bergenthum and Riischendorf (2007a)) and thus condition 3.3-(2) in Theorem 3.3 holds. 
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(c) (Smooth operator) Let A s be an operator as in Proposition 2.11, then, for transition operators of PII 
on C? c 2 (i/) resp. the condition 3.3-(l) reduces to show that T is a smooth operator. For instance, 

let f £ T° = T l~l VV 2 ( v ) and let the assumptions of Proposition 2.7 hold true. Then, for a smooth 
operator it holds that T a ^f £ C 2 (M d ). Further, due to Proposition 2.7, it holds that T Stt f £ C\(y). 
Since, df , d 2 f £ C\(y) we also have T Stt (df) and T s ^(d 2 f) £ E\{y). Hence, T S}t f £ VV 2 (^) C D(A S ). 


4 Applications 

4.1 Levy driven diffusion processes on Co(M d ) 

In the recent financial mathematics literature Levy driven diffusion processes have found considerable at¬ 
tention as stochastic interest rate models or stochastic volatility models. For a classical survey see Heston 
(1993). Recent results and developments in this field can be found in Poulsen et al. (2009) and the references 
therein. 

Let {L t )t >0 be an R -valued Levy process with local characteristics ( b, a, F). We consider the Levy driven 
diffusion processes defined as solutions of the following stochastic differential equation: 


dX t = $(. X t -,t)dL t 
Xq = x, x £ R d , 


where $ : R d x R + —»• R dxd is in C'^’ 1 (M <i x M+). 

In the present section we will make use of the (probabilistic) symbol. This opens a neat way to derive 
the structure of the generator of Levy driven diffusions and, more general, to the class of Feller evolution 
processes. 

Definition 4.1 Let X be a Markov process with right-continuous paths (a.s) and 

T s , t u{x) := E x ’ s u(X t ) = E(u(X t )\X S = x) 

for s < t. IfT s j is a strongly continuous ES on the Banach space Co of continuous functions vanishing at 
infinity, then X is called a Co-Feller evolution process. The process is called rich, if 

C c °° C T>(A S ) for every s > 0. 

Let us mention that in some parts of the literature Cf, is used instead of Co- This larger space is then 
endowed usually with the strict topology (cf. Van Casteren (2011), Section 2.1). Here, it is convenient 
to use Co as reference space. We will see in the proof of Lemma 4.6 below that in the setting (4.1) we 
obtain a Co-Feller evolution process for every Levy process. If we used Cf, instead, we would either loose 
the Banach space structure (using the strict topology which is defined via a family of semi-norms) or we 
would loose strong continuity of the semigroup in most cases (using the sup-norm). The following concept 
of a probabilistic symbol has proved to be useful in the time-homogeneous case already. We introduce this 
notion for the time-inhomogeneous case in order to derive path properties of the process. In the present 
article we will use it to derive the structure of the generator. 

Definition 4.2 Let X be an R d -valued Markov process, which is conservative and normal, that is P X (X 0 = 
x) = 1, and having right-continuous paths (a.s.). Fix a starting point x and a starting time t > 0 and define 
t = tO to be the first exit time from the ball of radius R > 0 after time t: 

t := t x r := inf {h > 0 : \\X x +h - x|| > R}. 
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The function p : R + x l d x l d —> C given by 


(4.2) 


pit , x, £) := — lim E x 

hiO 


e H x f+h~ x ’Z) _ i 


is called the (time-dependent) probabilistic symbol of the process, if the limit exists for every t > 0 and 
x,f £ R d independently of the choice of R > 0. 


Let us recall that in the above definition X[ +h = X min {t+ft,,r} denotes the process (which is started at 
time t in x) stopped at time r. 

At first sight it might be a bit surprising that it is possible to demand that the limit does not depend 
on the choice of R respectively r. Intuitively speaking, this is due to the fact that we are dealing with 
right-continuous processes. The symbol describes only the local dynamics of the process and this is not 
changed by using the stopping time. 

We need two auxiliary results in order to prove Theorem 4.5. The first one is the time-inhomogeneous 
version of Dynkin’s formula. It follows from the well-known fact that (M^ x ’*’“^>o given by 

t-\-h 

M [xa,u] :=u (x t+h ,t + h)-u(x,t)- j {A s + d s )u(X s ,s)ds 

t 

is a martingale for every u £ Pi t<s<t+h'D(A s ) ^ d?(d s ) with respect to every P x,t , x £ R d , t > 0. To our 
knowledge the most general version of this result can be found in (Casteren, 2011, Theorem 2.11). There it 
is formulated for the Cb case but it can be easily adapted to the Co case. For a function u which depends 
on the time-inhomogeneous process, but not on time itself, we obtain: 


Lemma 4.3 Let X be a Feller evolution process (Cb or Co) and r a stopping time. Then we have 

rA(t+h) 


E*’ 1 


A s u(X s ) ds = E x ’ t [u(X TA(t+ /,))] - u(x) 


(4.3) 


for allt > 0 and u £ Df< s <t+/t ^{A s ). 

In the previous lemma we have used that the class of martingales is stable under stopping. The following 
result has been established in Schilling and Schnurr (2010, Lemma 2.5). 


Lemma 4.4 Let K C R d be a compact set. Let \ '■ t R be a smooth truncation function, i.e., \ € 
C%° (M d ) with 

1 B 1 (o){y) < x(y) < is 2 (o) (y) 

for y £ R d . FuHhermore we define Xn(y) '■= x((y ~ x)/n ) and for (el 11 u x (y) := Xniy) 6 '’^’^ ■ Then for 
sufficiently large c > 0 we have for all z £ K 

|<P + y) - «a) - (y, V<p))i Bl(0 )(y)| < c • (1 a |y| 2 ) . 

Now we are in a position to show that the probabilistic symbol (which is easy to calculate) and the 
functional analytic symbol (which we will use below) coincide. 


Theorem 4.5 Let X = ( X t )t>o be a conservative cadlag Feller evolution process such that C^° C D(A S ) 
(s > 0). Then the generator A s \c°° is a pseudo-differential operator (cf. Section 5 below) with symbol 
—q(s, x, £), that is 

A.f{x) = -(2? T)~ d / 2 f e-^qfs, x, OHOdf (4.4) 

R d 

where the symbol q : R + x R d x R d —> C has the following propeHies for fixed s > 0: 
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• q(s, •, •) is locally bounded in x , 

• <z(s, •,£) i s measurable for every £ £ R d . 

• q(s, x, •) is a continuous negative definite function in the sense of Schoenberg for every x £ R d . 


T := t r := inf {/z > 0 : \\Xf +h - x\\ >R}. 
If x i-> q(s,x,£) is continuous, then we have 


lim E x,s 

hlO 


e i{X I +h ~X,£) 

h 


-q(s,x,f), 


(4.5) 


(4.6) 


independently of the choice of R respectively t, that is, the probabilistic symbol of the process exists and 
coincides with the symbol of the generator. 


Proof: As in the time homogeneous case, it is easily seen, that the operators A s fulfil the positive maximum 
principle. Therefore, they are pseudodifferential operators by Courrege (1966, Theorem 1.5). This is the 
functional analytic approach; now we complement this with the probabilistic approach. Let (x^)ngN be the 
sequence of cut-off functions of Lemma 4.4 and we write e^(x) := e l ( x & for x,f £ R d . By the bounded 
convergence theorem and formula (4.3) we see 


E x ’* - l) = Jim (E^[ x x n {Xf + MXf +h )e^{x) - l]) 

= e-i{x) lim E x ’* (xn( x [+h)ed x f +h ) - Xn( x ) e d x )) 


rA(H/i) 


= e-t(x) lim E x ’ t [ A s (x x e ( )(X s ) ds 

n—> oo J 

t 

rA(t+/i) 

= e-t(x) lim E x ’ t [ A s (xf l e i )(X s _) ds. 

n—> oo J 


The last equality follows since a cadlag process has a.s. a countable number of jumps and since we are 
integrating with respect to Lebesgue measure. Using our Lemma 4.4 and writing the operator A s in integro- 
differential form, we obtain for all z £ Br{x) 


A s (Xnei)(z) 

< c x{\b(s,z)\ + i E k*mi + / (iAtei 2 )n»^.<i»))(i+i{i 2 ) 

»*> 

<c' x sup sup |p(s,z,g)| 
z£Br(x) l€l<! 


where c x and c' x are positive constants only depending on %. The last estimate follows with techniques 
described in the appendix of Schilling and Schnurr (2010). The function p(t,x,tf) is locally bounded since 
it is continuous. By definition of the stopping time r we know that for all t < s < r A (t + h) we have 
X s _ £ Br{x). Therefore, the integrand A s (%“e{)(X s _), t, < s < t A (t + h) appearing in the above integral 
is bounded and we may use the dominated convergence theorem again. This yields 

E - l) = e_^(x)E x,t f lim A s ( X x n e^(z)\ z=Xe _ ds 

\ J J n—> oo 

t 
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tA(£+/i) 

= -e- i {x)E x ' t J effz)p(s,z,f)\ z = Xs - ds. 

t 


The second equality follows from classical results due to Courrege (1966, Sections 3.3 and 3.4). Therefore, 

E*,t ( e K x t+h~ x <i) — i) 

lim--- 

hiO h 


= —e-e(x) limC 11 ’* 
^ ft.4.0 


= —e-e(x) lim E x,t 
^ h.4.0 



since we are integrating with respect to Lebesgue measure. The process X T is bounded on the stochastic 
interval [[0, r[[ and (s, x) >—> p(s, x, f) is continuous for every f £ R d . Thus, again by dominated convergence, 

E x ’ t ( e i ^ x I+h~ x ’i) — 1 ) 

lim--- - -= -e-t(x)es(x)p(t,x,£) = -p{t,x,g). □ 

/40 u 


In fact, above we have been more general then needed in the context of equation (4.1). For bounded 
coefficients (and hence a bounded symbol) it is not necessary to introduce a stopping time. In this case the 
limit without stopping time, that is, 


— lim E x,t 

/i.j.0 


J(X t +h-x,0 _ 1 
t 


would have been sufficient and coincides with p(t,: r,£) and q(t,x,£) above. For future reference, however, 
we have included the general case. Now we use the above result on Levy driven diffusions. 


Lemma 4.6 The unique strong solution of the SDE (4.1) Xf(u>) has the symbol q : M+ x R d x R d C 
given by 

q(t, x, £) = p(t, x, 0 = ip($ T (x, 

where is the coefficient of the SDE and tp the symbol of the driving Levy process. Hence, the (family of) 
generators on T>(A S ) D C'^R 6 *) can be written as 


A sf( x ) = l Y, M(^s)) J ’ fc ^^-( a; ) + ^(6$(x, S )) J ^ : (a ; ) 


j,k=i 


l=i 


+ 


J (f{x + ${x,s)y) - f{x) ~Y §~^ x )(^(x, s)xc(y)Y^F(dy), 

i=i 3 


(4.7) 


where F is the Levy measure corresponding to L which integrates (1 A \y \ 2 ) and \c is a cut-off function. 

Proof: It is a well-known fact that the unique solution of a Levy driven SDE with Lipschitz coefficient is a 
Markov process (cf. Protter (1977, Theorem 5.9). In order to carry the solution, the stochastic basis on which 
the Levy process is defined has to be enlarged in a canonical way (cf. e.g. Protter (2005, Section V.6)). Next, 
one has to show that the solution is a universal time-inhomogeneous Markov process, that is, the transition 
function does not depend on the starting point. This can be done similarly to Schnurr (2009, Theorem 2.47). 
The Co-Feller property follows as in Schnurr (2009, Theorem 2.49). 
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The calculation of the symbol works as in the proof of Schilling and Schnurr (2010, Theorem 3.1). We 
sketch the proof here in the case d = n = 1 in order to emphasize analogies and differences: let r be the 
stopping time given by (4.5). Fix x, £ e M. We apply Ito’s formula for general semimartingales to the 
function e^(- — x) = exp(f(- — x)£): 


- E x,t ( e i ( x * T-a: )£ - l) 

t~\~h 

= ^-E x ’ t (^J i£e i{x Z-- x)( dX J-i J fe KX --~ x)( d[X T ,X T ] c s 

t+ t+ 

+ e“ ix€ ( eiX ‘ ? - eiX °- 5 - e**'" 5 AX s) ) ■ 


t-\-h 


(4.8) 


For the first term we get 

t-\-h 

j E x,t J (i£e i(x Z-~ x) ^ <LYJ 

t+ 


= - E x ’ 
h 

t+h 






7 < 

'i£ e i( - x *- 

-«) d / 


r “) 1 [[0,- 

r]](-,r)dZ, 

t+ 


Vo 




= - E x ’ 
h 

t+h 






7 < 

i£ e^ x “- 

_a:)5 3>(X s _ 

> s -)![[o, 

,r]] (■»»)) 

dZ s 

t+ 






= - E x ’ 
h 

t+h 






7 ' 



> s -)l[[o : 

,d] (■>*)) 

d(bs ) 


(4.9) 


t+ 

t-\-h / 

+ t E x ’ t f (^e i ( x J--^«$(X s _, S -)l [[0iT]] (., S ))d AZ r l { | AZr |> 1} 

t+ \0 <r<s 

where we have used the well-known Levy-Ito decomposition of the driving process. Since the integrand is 
bounded, the martingale terms of the Levy process yield again martingales whose expected value is zero. 
For the drift part (4.9) we obtain 


t~\~h 

l EX,t J (^•e i ^-- x )«4>(X s _, S -)l [[0jT[[ (., S )6) ds i£M>(x,t), 

t+ 


\ EX,t j (^■e < W--*> € $(X a _,a-)l [[0lT[[ (-,a)6) ds 

t+ 

t+i 

= ■ E^ j ( e «-*)€$(x afc) sh) 1 [[0) t[[(-, shj) ds 

t 

The other parts work alike. In the end we obtain 

p(t,x,£) = -ib($(x,t)€) + l($(x,t)£)a($(x,t)£) 

y^o 
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M)0- 


(4.10) 


Note that in the multi-dimensional case the matrix has to be transposed, that is, the symbol of the 

solution is ip(Q T (x, f)£). We already know that the probabilistic and the classical symbol coincide. Hence, 
some Fourier analysis gives us the structure of the generator. □ 


Remark 4.1 We need the boundedness of $ in order to derive the Feller property. The calculation of the 
symbol works even for unbounded coefficients (e.g. Lipschitz continuous <f>). 

Now we can prove the comparison result for Levy driven diffusion processes. We consider T to be a subset 
of Co(R d ) and (L^) t >o, i = 1,2 to be a Levy process with local characteristics (b^\ a^\ F^). Further, 
let X and Y be Levy driven diffusions defined by the stochastic differential equation (4.1) associated to 
(dh 1 ), L^ 1 )) and (& 2 \ L^) with corresponding infinitesimal generator A s and B s respectively. Denote their 
evolution systems by S = {S s ,t)t> o and T = (T s>t ) t >o, respectively. 

Corollary 4.7 (Comparison of time-inhomogeneous Levy driven diffusion processes) 

Assume that Xo = Yq and 

(1) / € T implies s H > T s ^f is right-differentiable for 0 < s < t, 

(2) X is stochastically monotone w.r.t. T i.e. f € T implies S s jf G T for all s < t, and 

(3) A s f < B s f [P A '°] for all f G T and s <t. (4.11) 

Then 

X t <?Y t , t> 0. (4.12) 

Proof: We have seen above already that for every choice of and every Levy process L^) the solution of 
(4.1) is a Co-Feller evolution process. Since (Co, ||-|| ) is a Banach space we can directly use Theorem 3.3 
and obtain S s jf < T s jf for all / G J 7 , s < t. By using Xo = Yq this yields for / £ J and t > 0 

Ef(X t ) = E 0 [E[f(X t )\X 0 = x]] 

= E 0 S 0 , t f(x) 

< E 0 T 0tt f(x) 

= Eo[E[f(Y t )\Y 0 =x]]=Ef(Y t ), 

where Eq is the expectation with respect to P x °. □ 

Corollary 4.7 shows that ordering results for Levy driven diffusion processes for JcC 0 2 are implied by a 
stochastic monotonicity assumption of one process and comparability of the local characteristics. In Section 
4.3 we consider Levy driven diffusion processes as inhomogeneous Markov processes on C%(v) in order to 
derive comparison results w.r.t. unbounded function classes. Some further applications of these ordering 
results derived in a similar way are listed in the following remark. 

4.2 Processes with independent increments 

Let (L t )t> o be a PII such that E\L t \ 2 < oo for all t G M+. Its characteristic function is given by Ee l ^ lLt ' 1 = 
exp (f* 9 u (if)du) with cumulant function given in (2.18). We consider the case of PII with Levy measures 
Fs l \i = 1,2 which integrates (|y| A \y\ 2 )- We choose the cut-off function Xc as the identity. Denote for s > 0 
the local characteristics of L^> by (b s . cr^'K Fg'^). The corresponding infinitesimal generators are A s and B s 
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respectively. If the assumptions of Proposition 2.7 are fulfilled the family (T s>t ) s < t of operators defined in 
(2.11) is a strongly continuous ES on C\(v). The corresponding infinitesimal generator (H s ) s >o is given in 
(2.19) defined on D(A S ). For this instance we state a comparison result for the convex order < cx which is 
generated by 

j? x = 7 a ns b nw 2 M, 

where B b is the function class of 6-bounded function and b is chosen as in Remark 2.8. Similar results hold 

for :<]<:>:■ ^smi ^st; A ism aS well. 


Corollary 4.8 (Comparison of PII w.r.t. -F cx ) 

Assume L ^ = Ll?' 1 for some r < t. If ( S s> t ) and (T s y) are smooth operators and if (S St t) leaves Bi, 
invariant, then the comparison of the local characteristics 

< P sd <?i 2) (4.13) 

J f(x)F^ 1) (dx) < J f{x)Ff>{dx) (4.14) 

R d R d 

for all s < t and all f £ F cx with /(0) = 0 such that the integrals exist, implies that 

<cx l[ 2 \ t > 0. (4.15) 

Proof: Let / £ Since (S St t) leaves By, invariant we obtain by Remark 3.4-(c) and the smoothness of 
(S S 't) that S s jf £ VV 2 (^). The invariance in F cx follows by Remark 3.4-(b). Similarily the smoothness of T Sjt 
and Proposition 2.11 ensure that T Sjt / £ V(B S ) for all s > 0 (see Remark 3.4-(c)). We have choosen above 
the cut-off function xc as the identity, since the Levy measures fulfil / K d(|j/| A \y\ 2 )Fs l \dy) < oo, i = 1,2. 
Thus the jump part has the form 


J f(x + y)- f(x) £ y-(x){y)WFV{dy) 

R d j=1 3 

= J f(x + y) ~ f(x) -J2^-{x)(y) b) l {lyl<1} F^\dy) 

R<i ■ 7 = 1 3 

- j ^Z§[:(x)(y) ij) t { \ y \>i}F^(dy). 


From this representation we note that both terms belong to C\(y) which can be seen by using similar 
arguments as in the proof of Lemma 2.10. Thus, the infinitesimal generator has the form as in (2.19). 
Hence, for / £ Ff? x and s < t we obtain 


(■ B s - A s )f(x) 



j,k=l 


d 2 f 

dxjdxk 




+ J (f(x + y)~ f{x) 

R d 


§^r.-i x )i y ) 3 ) ( F s 2) - ^i 1} ) (dy)- 

4=1 3 


(4.16) 


(2\j,k 

Due to the positive semidefiniteness of (a s ; —a\ )j,k<d for fixed s > 0, its spectral decomposition is 
given by (JT <d A iejef ), where the eigenvalues A* are non-negative and e = (ej,..., ef) denote the eigenvec¬ 
tors. Convexity of / implies 


d 


d 


d 2 f 


2 dxjdxk 

i—l j,k—l J 


{xy^i > o. 
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Moreover, the function y H> f(x + y) — f(x) — Y^j=i [x)y^ is convex, since it is a sum of a convex function 
and a linear function. Now using the comparison of the Levy measures in (4.14) we get ( B s — A s )f(x) > 0. 
Thus, the assumptions of Theorem 3.3 are satisfied and we obtain 

EfiL^) = E r E{f{L^)\L^ =x)< E r E{f{L^)\L^ = x) = Ef{Lf\ 


where E r is the expectation with respect to P 


4 1 ’ 


□ 


Remark 4.9 The following gives for some function classes comparison conditions of the local characteristics 
which imply stochastic ordering of the Levy driven diffusion processes in the case of Levy measures, F ^>, i = 
1,2 which integrate (|y| A \y \ 2 ). For d x d matrices A , B with real entries, the positive semidefinite order 
A < P sd B is defined by x T (B — A)x > 0 for all x £ R d , where x T is the transpose of x. Here, we make use 
of the explicit representation of the generator given in (4.7). 

If : R d x R+ ->• R^ xd , i = 1,2 and $b) < ps d then sufficient conditions for (4.11) are: 


Ordering 

Drift 

Diffusion 

Jump 

Tst 

6« < fe (2) 

crb) — (jfA 

C? 

fe. 

II 

X sm 

6b) = 6 (2) 

fjb) < d 

F W <sm f (2) 

T 

J ism 

< b m 

fjb) < d (jO 

1 ) < /?( 2 ) 
—sm 1 


where crb) <d cr ( - 2 - ) means pointwise ordering cr*- 1 ) < cr^ 2 \ — a-( 2 ) JJ , j < d. For every generating 

function class, the latter fact can be justified as follows: Let, for instance, f £ F st . Then, the first or¬ 
der derivatives of f are non-negative. Moreover, let &b) < b^ 2 ' 1 <jb) = o^ 2 ) and F^ = F^ 2 \ hence the 
corresponding generators in (4.7) are ordered, too. Then, the comparison result follows from Corollary j.10. 

4.3 Levy driven diffusion processes on 

We reconsider Levy driven diffusion processes as inhomogeneous Markov processes. In comparison to Corol¬ 
lary 4.7 where these processes were considered as strongly continuous evolution systems on Co(R d ) we now 
consider them as strongly continuous evolution systems on (y) assuming for both processes the kernel 
assumption ( K) as well as (C) and (£>). We state a comparison result for Levy driven diffusion w.r.t. <jro, 
where T is from (1.3) - (1.7) and 

J° = J n 6 b n W 2 (v). 

Now let (L^) t >o, * = 1,2 be a Levy process such that E\L t \ 2 < oo, t > 0, with local characteristics 
(2>W, o -(0 ) fW) such that F^\ i = 1,2 integrates (|y| A \y\ 2 ). Further, let X and Y be Levy driven diffusions 
defined by the stochastic differential equation (4.1) associated to (f&b), L^) and ($( 2 \ 4 2 )) respectively. 
Denote their infinitesimal generator by A s and B s respectively defined in (4.7) and their evolution systems 
on L\{v) by S = (5 s>t ) t >o and T = (T Sit )t>o- 

Extending the ES of X and Y onto C 2 {v) by Proposition 2.9 it is clear that the form and the mapping 
behaviour of the associated infinitesimal generators (see equation (4.7)) for f £ VV 2 ffi) is deduced similarly 
as in Lemma 2.10 and Proposition 2.11. 

Corollary 4.10 (Comparison of Levy driven diffusion processes) 

Assume that X and Y possess the transition kernels (P f x ) t >o and (P^ ) t > o, which satisfy ( K ), ( C ) and ( D ). 
If X 0 = Y 0 and 

(1) / £ T implies s i-> T S)t f is right-differentiable for 0 < s < t, 

(2) X is stochastically monotone w.r.t. F°, and 
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(3) A s f < B s f a.s. for all f G T and s < t, (4-17) 

then 

X t <? Y t , t > 0. (4.18) 


5 Pseudo-differential operators on Sobolev Slobodeckii spaces 


Finally, we make some remarks on a method to compare general Markov processes based on Sobolev Slo¬ 
bodeckii spaces and the theory of pseudo-differential operators. We relate this method to strongly continuous 
evolution systems as used in this paper. The main reference for this final section is Bottcher (2008) who gave 
a representation result for pseudo-differential operators on Sobolev Slobodeckii spaces similar as in Theorem 
3.1. His result is based on developments on pseudo-differential operators in Eskin (1981), Hoh (1998) and 
Jacob (2001, 2002). 

For reM the Sobolev-Slobodeckii space is defined as 

n r (R d ) := {u G «S'(R d )|u G Aoc(R d ) s-t. |M| r < oo} (5.1) 


with 


Nl?= / |fi(0l 2 (i + l£l) 2r dC. 

R d 


(5.2) 


where <S , (R d ) is the dual space of the Schwartz space <S(R d ) and u denotes the Fourier transform of u. The 
space 77 r (R. d ) endowed with the scalar product (•, •) r , 


(«, v) r := J «K)t)(0(l + kl ) 2r df, u, v G ?T(R d ) 


is a separable Hilbert space (see Eskin (1981, Theorem 4.1)). Moreover for r > 0 ( / H r (M d ), (•, -) r ) is isomorph 
to 


{u G £ 2 (R d ,C)| J |u(0| 2 (! + l^l) 2r ^ < oo} 


endowed with the same scalar product. The space W( R d ) is an interesting example of a Banach function 
class which uses regularity conditions on the elements but allows unbounded functions. 

Elements u of 'H r (M d ) fulfil certain integrability and regularity assumptions since the existence of higher 
moments of u in (5.2) corresponds to higher order regularity of the function. For instance, u G A)!!) is 
continuously differentiable if the first moment of u is integrable, that is if J |£||u(£)|d£ < oo, and in this case 
we have 





A pseudo-differential operator q(t, x, D) is defined on a suitable space by 


q(t,x,D)f(x) := (2tt) d/2 f e 

M d 


(5.3) 


if the right-hand side exists. The function q which defines the operator is called the symbol of the pseudo¬ 
differential operator. 

In Bottcher (2008, Theorem 2.3 ) a representation result for solutions of evolution equations on W{ R d ) 
is established. For an application of this result it has to be checked that the transition functions T Sjt and the 
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infinitesimal generator A s (corresponding to a time-inhomogeneous Markov process) are pseudo-differential 
operators with symbols in suitable ’symbol classes’, i.e. the symbols have sufficient regularity and growth 
behaviour. Under this condition on the symbol it is verified in Theorem 1.8 in Bottcher (2008) that the 
transition function T s>t maps "H' (R d ) into itself. This mapping property is based on general theory of 
pseudo-differential operators as described in Eskin (1981), Hoh (1998) and Jacob (2001, 2002). 

Our representation result in Theorem 3.1 can be applied to this context if the transition operator T Sjt is 
a strongly continuous ES on H r (M. d ), that is, it maps "H r (R d ) into "H r (R d ). In other words, it is sufficient to 
prove an analogous result to Proposition 2.4 for W(\R d ). 

We can apply our representation Theorem 3.1 and get as corollary a result related to Bottchers repre¬ 
sentation result without any growth and regularity conditions on the symbol of the infinitesimal generator. 

Corollary 5.1 Under the assumptions of Theorem 3.1, let ( T Stt )s<t he a strongly continuous ES on r H r (M d ), 
then the solutions Ft of the weak evolution problem (3.2) on TL r {M. d ) have the representation 

t 

F t (s) = T', t F t {t ) - J T Str G(r )dr. (5.4) 

S 

For the application of Corollary 5.1 first note that any transition operator (T s>i ) of a Markov process 
(X t ) t > o has a representation as a pseudo-differential operator. 

For the proof let / G % r (M d ), r > 0 and consider a time-inhomogeneous Markov process X = (X t ) t >o 
with transition function ( P s ,t)s<t ■ Denote by F(f) the Fourier transform of / and by F^ 1 the inverse Fourier 
transform. Then 

T s , t f(x) = E s , x (F- 1 (F(f))(X t )) 

=E Aw^I e ~' <x, ' i>Fmod( ) 

R d 

= ( 27r )rf/2 J e ~ i(x,i)< ls(t,x,Of(Od^ (5.5) 

R d 

where q s (t,a :, £) := E(e~ l( ' Xt ~ x ’^\X s = X s ). Thus (T S)t ) is a pseudo-differential operator with symbol 

Qs(t,x,£). 

For applications of Corollary 5.1 to comparison results we have to establish that T s>t : W(R d ) —» W(R d ). 
This can be done in the case of PII, which are particularly suitable for the pseudo-differential operator 
approach. 

Example 5.1 (Transition functions of PII on H r (M. d )) 

Let X = (Xt)t> o he PII with characteristic function £ , (e _ ^ x *~ x,s ’^) = exp (f*6 u (—if)du) =: /x s ,*(—£)■ 
Then from the representation in equation (5.5) we obtain for f G TL r { R d ), that 

Ts,tf(x) = J (5-6) 

R d 


Hence the characteristic function fi s ,t{—£.) is the symbol of T sit . It is not difficult to see that the transition 
operator of a PII is a pseudo-differential operator with bounded C°° -symbol, if all absolute moments of the 
Levy measure exists. Thus the symbol lies in the symbol class Sq. In consequence the conditions on the 
symbol for the mapping property T Stt : 'H r (M d ) H r (R. d ) in Theorem 1.8 of Bottcher (2008) are fulfilled 
and Corollary 5.1 gives the necessary representation result, result Corollary f.8. 

As consequence of Corollary 5.1 the comparison theorem (Theorem 3.3) allows us to state a comparison 
result for PII for function classes T C 7t r (R d ). Note however that by this approach, based on H r (M. d ), we 
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need the strong condition of existence of all absolute moments of the Levy measure, while our comparison 
result based on modified C\{v) spaces as in Corollary 4-8 is established under much weaker conditions. 

For general Markov processes with symbol q s (t, x,f) dependent on x it is not easy to check that its symbol 
has suitable regularity and growth behaviour. Hence this method is difficult to apply for time-inhomogeneous 
Markov processes as for example for time-inhomogeneous diffusions. Thus it seems that the approach in this 
paper based on our representation result for strongly continuous evolution system in general Banach spaces 
is more flexible and easier to apply in examples. 

Acknowledgement: AS would like to acknowledge financial support by the DFG (German Science Foun¬ 
dation) for the project SCHN 1231/2-1. 
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